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Abstract 

We consider semi-stable, radially symmetric, and decreasing solutions 
of —A p u = g(u) in the unit ball of M n , where p > 1, A p is the p-Laplace 
operator, and g is a locally Lipschitz function. For this class of radial so- 
lutions, which includes local minimizers, we establish pointwise, L q , and 
W 1,q estimates which are optimal and do not depend on the specific non- 
linearity g. Among other results, we prove that every radially decreasing 
and semi-stable solution u belonging to W 1,P {B\) is bounded whenever 
n < p + Ap/ (p — 1). 

Under standard assumptions on the nonlinearity g(u) = Xf(u), where 
A > is a parameter, it is proved that the corresponding extremal solution 
u* is semi-stable, and hence, it enjoys the regularity stated in our main 
result. 
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1 Introduction 



This article is concerned with react ion- diffusion equations —A p u = g(u) involving 
the p-Laplace operator. We consider a class of radially symmetric solutions: those 
solutions which are decreasing and semi-stable. This type of solutions include 
local minimizers, minimal solutions, extremal solutions, and also certain solutions 
found between a sub and a supersolution. We establish sharp pointwise, L q , and 
W 1,q estimates for solutions in this class. In particular we show that semi-stable 
radial solutions enjoy better regularity properties than general radial solutions. 
In addition, our results do not depend on the specific form of the nonlinearity in 
the reaction term. More precisely, our pointwise and L q estimates hold for every 
locally Lipschitz nonlinearity g, while the W l,q estimates hold for nonnegative g. 
All the results obtained in this paper were obtained by the first two authors in 
[5] for the Laplace operator, namely p — 2. 

We consider radially symmetric and decreasing solutions u £ W 1,p (Bi) of 

- div(\Vu\ p - 2 Vu) =g{u) in B x \ {0}, (1.1) 

where p > 1, B\ is the unit ball of K n , and g : R — > R is a locally Lipschitz 
function. By a radially decreasing function u we mean a function u such that 
u = u(r) and u r (r) = (du/dr)(r) < for all r £ (0,1), where r = \x\ and u r 
denotes the radial derivative. 

We do not assume u(l) = or any other boundary condition. Nevertheless, 
since u £ W 1,P (B\) is radial, we have u £ L^ C (B i\{0}) by the Sobolev embedding 
in one dimension. Hence, using known regularity results for degenerate elliptic 
equations (see [2"Tj). we have that in fact u £ C^(Bi \ {0}) for some (3 £ 
(0, 1). However, u may be unbounded at the origin. Our main result establishes 
estimates for semi-stable solutions in the whole ball B\. Thus, this may be 
regarded ClS db result on removable singularities. 

Consider the energy functional 

E a (u):=- I \Vu\ p dx- [ G(u)dx, (1.2) 
V Jn Jq 

where G' = g and Q is a smooth bounded domain of R n . We say that a radially 
decreasing function u £ W 1 ' P (B 1 ) is a radial local minimizer of fll.2p with Q = B t 
if for every 5 > there exists eg > such that 

for all radial functions £ £ Cl(Bi \ B$) (i.e., with compact support in B\ \ B$) 
satisfying ||£||ci — £s- Note that every radial local minimizer u is a solution 
of (11.11) . Moreover, u is semi-stable in the sense that the second variation of 
the energy at u is nonnegative. The following definition makes this precise — see 
Remark 1 1.71 for more comments on this, as well as [7j for the corresponding setting 
in the nonradial case. 
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Definition 1.1. Let u G W 1,p (Bi) be a radially symmetric solution in B\ \ {0} 
of (11. ip such that u r (r) < for all r G (0, 1). We say that u is semi-stable if and 
only if the second variation of energy Q at u satisfies 

Q{Z)-= I {(p-l)kr 2 |^| 2 -9'Wdrfx>0, (1.3) 

for every radially symmetric function £ G C^(I?i \ {0}). 

As we will see later, the class of semi-stable solutions includes not only local 
minimizers but also the minimal and extremal solutions of the classical problem 
(1.12a iP ) below, which motivated our work. 

To state our estimates, we define exponents qk, for k = 0,1, by 



1 1 2/72-1 k - 1 2 4p 

— for n > p + 



q k p np\ p — 1 n np p — 1 (14) 

q, := +oo tor n < p H . 

p-1 

One can easily check that p < q k < +oo in all cases. In addition, 



ti — p in — 1 A.p 

> 1 + \ if and only if n > p H . (1.5) 

2 \ P — 1 p — 1 

Hence, both exponents qo and q\ are well defined. Moreover, q and qi are finite 
if n > p + Ap/ (p — 1). 

Our main theorem states sharp regularity results and pointwise estimates 
for every semi-stable radially decreasing solution u of (11. II) . It also establishes 
estimates for its radial derivative u r . Our result extends Theorem 1.5 and part of 
Theorem 1.8 of [5] (in which p — 2) to the general case p > 1. In particular, if we 
set p = 2 in the following theorem, one recovers the results of [5] for the Laplace 
operator. The statements in the following theorem are optimal in several respects 
discussed below, and are still open in general domains even for p = 2. Indeed, in 
the nonradial case, they are known to hold only for g(u) = Xe u , g(u) = A(l + u) m 
with m > p — 1, and nonlinearities g close to them in certain senses — see the 
comments following Theorem 11.31 for more details. 

Theorem 1.2. Let g be a locally Lip schitz function andu G W 1,p (Bi) be a semi- 
stable radial solution in B x \ {0} of (11.11) satisfying u r {r) < for all r G (0, 1). 
Then: 

(a) If n < p + 4p/ (p — 1) then u G L°°(Bi). Moreover, 
where C n ^ p is a constant depending only on n and p. 
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(b) If n = p + 4p/ (p — 1) then u G L q (B\) for all q < +00. Moreover, 

\u{r)\ < Cp||w|| w i, P ( Bl) (|logr| + 1) in Si, (1.6) 

where C p is a constant depending only on p. 

(c) If n > p + Ap/ (p — 1) and q < q , then u G L q (Bi) and 

IM|l<7(Bi) < Cn,p,g||'u||wl.P(_B 1 ), 

where C n ^ PA is a constant depending only on n, p, and q. Moreover, 



.I(n-2,/SEf-p-2 



|«(r)| < C niP ||n|| wl , P(Bl) r p V VH y(|logr|p +1) mSx, (1.7) 

where C nyP is a constant depending only on n and p. 
(d) Assume that g is nonnegative. Then: 
(dl) W^e /iove 



l|v«|U P(Bl) < c^|||(ti-ti(i))^ l ||^ fli) + \\g(u)\\i^ Bi) 

for some constant C niP depending only on n and p. 
(d2) u G W 1,q (Bi) for every q < g 1; and 

\\u\\ w ^(Bi) < G ifq<qi, (1.9) 

where C is a constant depending only on n, p, q, and on upper bounds for 

IMUh-Bi) and 9- 

(d3) If n >p + 4p/ (p — 1) t/ien 



'li 



|t*r(r)| < C ni p||w|| w i, P ( Bl) r P V VP 1 7|l gr|p in B 1/4 , (1.10) 
where C n ^ p is a constant depending only on n and p. 

Examples 11.41 and 11.51 below on the exponential and power nonlinearities show 
the sharpness of the previous regularity results. Indeed, the functions 

V 

u(r) = — plogr and u(r) —r m-O- 1 ) — 1 (1-H) 

are solutions of (11.11) for, respectively, g(u) = X*e u and g(u) = A*(l + u) m with 
certain values of A* = X*(n,p, m). Using weighted Hardy inequalities, one can 
easily find the ranges of values of n, p, and m for which these solutions are 
W l,p (Bi) semi-stable solutions (see [7] and Examples 11.41 and 11.51 b elow for more 
details). The ranges obtained in this way for the exponential nonlinearity show 
the sharpness of condition n < p + 4p/(p — 1) in Theorem 11.2( a). as well as the 
optimality of the pointwise estimate (jl.6p in part (b). The ranges obtained for 
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the power nonlinearities give the optimality of the exponents go and q\ in parts 
(c) and (d), as well as the sharpness of the pointwise bounds (11.71) and (11.101) 
except for the factors | logr| 1//p in them. 

The proof of Theorem 11.21 was inspired in the proof of Simons theorem on 
the nonexistence of singular minimal cones in W 1 for n < 7 (see [6] for further 
details). The main point of the proof is to obtain the following key estimate 



(1 

\v \ p r~ 2a dr< — ^ HV?dl p 

1 r| ~ (n - 1) - (a - lf{p - 1) 11 VM| W0 




— from which most of our results follow — for those exponents a making the 
denominator in the above right-hand side positive (see Lemma [2T3T) . This estimate 
is established by taking £ = u r r] as a test function on the semi-stability condition 
(II. 3p . being 7] essentially a negative power of r. Lemma [2721 shows that, with this 
choice of £ in (jl.3p . the term g'(u) in (II. 3p disappears. This is the reason why 
our main estimates do not depend on g. 

As an application of Theorem 11.21 and the ideas behind its proof, we consider 
the following problem 

= Xf(u) in Bi, 

> in B u (1.12 A , P ) 

= onaBi, 

where A > 0, and / is an increasing C 1 function with /(0) > and 

li m ZW =+00 . (1.13) 

Problem (1.12 Ap ) is studied in [7] for general smooth bounded domains Q 
of M n . It is proved there that there exists a parameter A* 6 (0, oo) such that 
if < A < A* then (1.12\ tP ) admits a minimal solution u\ G C 1 (fi). Here 
minimal means smaller than any other supersolution of the problem. Moreover, 
for A > A* problem (1.12a iP ) admits no regular solution. We also have that 
for every < A < A* the minimal solution u\ is semi-stable (in the sense of 
Definition 11.11 when Q = B\). Consider the increasing limit 

u* := lim w A . (1.14) 

ATA* 

In contrast with the case p = 2 involving the Laplacian, it is not always clear that 
the limit u* is a weak solution of (1.12 AiP ) for A = A*. When one can establish 
that u* is a weak solution (this may depend on the assumptions on n, p, Q, and 
/), it is called the extremal solution. 

In our next result we prove that, for Q = Bi, the limit u* is actually a semi- 
stable radially decreasing energy solution of (1.12 A * p ). As a consequence, it enjoys 
the same regularity properties as the ones stated in Theorem 11.21 In particular, 
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we obtain that the extremal solution u* is bounded if n < p + 4p/(p — 1). As 
we will see below, the extremal solution is unbounded if n > p + 4p/(p — 1) 
when f{u) = e u . Therefore, this range of dimensions is optimal. The optimal 
dimension ensuring the boundedness of the extremal solution in general domains 
remains unknown for general nonlinearities / under the assumptions above (and 
even with the additional hypothesis that / is convex). 

Theorem 1.3. Let f be a positive and increasing C 1 function in [0, +oo) satis- 
fying (11 . 131) . For A G (0, A*) let u\ be the minimal solution of (1.12^), and let 
u* be defined by ( TTTT41) . 

Then, 

\\u\\\w^(Bi) + \\f(ux)\\ L i( Bl ) < C (1.15) 

for all A G (0, A*) and some constant C independent of A. Moreover, u* G 
W l,p (Bi) andu* is a semi-stable radially decreasing energy solution of (l.l2\* tP ). 
As a consequence, u* has the regularity stated in Theorem \\.2\ In particular, 

u* G L°°(Bi) if n<p+^—. (1.16) 

p — 1 

As mentioned in the statement of the theorem, for problem (1.12\ }P ) we con- 
sider energy solutions. That is, nonnegative functions u G W ,p (Bi) such that 
f(u) G L 1 (B 1 ) and 

/ \\7u\ p - 2 \7u- Vtpdx = [ Xf{u)cpdx for all (p G Cl(B%). (1.17) 

Note that, by a standard density argument, (11.171) holds for every ip G W^' p (Bi). 

The literature studying minimal and extremal solutions in general domains of 
M. n is extensive. Crandall and Rabinowitz [8] and Mignot and Puel [22] considered 
the case of the exponential and power nonlinearities in general domains for p = 2. 
They proved that u* is an energy solution for every dimension and that it is 
bounded in some range of dimensions. These results in general domains were 
extended for every p > 1 by Garcia-Azorero, Peral, and Puel [TTl [18] for the 
exponential case, and by Ferrero [16] and Cabre and Sanchon |7J for power type 
nonlinearities. 

An optimal regularity result for the extremal solution is not known for general 
nonlinearities and general domains, even for p = 2, with the exception of the 
radial case Q — B\. In [5J, Cabre and Capella obtained the optimal regularity of 
the extremal solution for every locally Lipschitz nonlinearity / when the domain 
is a ball and p = 2. In particular, they proved that the extremal solution is 
bounded if n < 9. In this paper, we obtain an analogue optimal radial result for 
all p > 1. 

In the nonradial case, [23] and [2] contain the best results for p = 2. In [23J, 
Nedev proves for p = 2 that the extremal solution is bounded if n < 3 — i.e., 
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n < 4 in (j!.18p — whenever / is an increasing and convex function such that 
/(0) > and f II . 1 3 j) holds. The same result up to dimension n < 4 is proved by 
the first author in [2j for p = 2 and Q convex, without the convexity assumption 
on /. On the other hand, the third author extends in [27] the work of Nedev [23] 
to the case p > 2 and establishes 

u*eL°°(Q) if n<p+^—. (1.18) 

p — 1 

Note that for p = 2, we recover Nedev' s condition n < 4. Note also the existing 
gap between the dimension in (11.181) for general Q and the ones in our result for 
Q — B% — note the presence of the factor 4 in ( 11. 161) in contrast with (11.181) . 

We refer to [101 E] for surveys on minimal and extremal solutions and to 
[H HH [121 H31 HH HSl H91 EDI I2S] for other interesting results in the topic of 
extremal solutions — also for problems involving other operators, such as the 
bilaplacian, or fractional Laplacians related to boundary reactions. 

The key point in the proof of Theorem II .31 is to establish the uniform estimate 
(11. 1 5 [) on the W 1,p -norm and the L 1 (S 1 )-norm of u\ and f(u\), respectively. This 
is accomplished using the "superlinearity" hypothesis (11.131) on /, together with 
the radially decreasing character of the solutions. The W 1,p bound is then used 
to show that the limit u* is an energy solution of (1.12a* iP ). Using again that the 
u\ are semi-stable, we prove that u* is semi-stable. 

In [7] , the semi-stability property of the extremal solution u* was proved in 
the case p > 2 for general domains in M. n and nonlinearities having its growth 
comparable to a power. This property was obtained taking the limit as A —>■ A* in 
the semi-stability condition for the minimal solutions u\. In the case 1 < p < 2 
this argument does not apply for general domains, since we have no control on 
the set where the gradient Vw* vanishes. However, in the radial case this set is 
the origin, and thus we obtain here the result of [7J also for 1 < p < 2. 

Next, we introduce two explicit examples which show the optimality of The- 
orems 11.21 and 11.31 

Example 1.4. In [T7] [T8] Garcia- Azorero, Peral, and Puel considered problem 
(1.12a iP ) for f(u) = e u and a general bounded domain Q. They proved that if 
n < p + 4p/(p — 1) then the extremal solution u* is bounded. On the other hand, 
if n > p + 4p/ (p — 1) and Q = B\ they show that 

u*{r) = — plogr and A* = p p ^ l {n — p), (1-19) 

and hence, that the extremal solution is unbounded. This example shows that 
the conditions on the dimension n in terms of p in Theorem II .2( a). (b) and the 
pointwise estimate (I1.6P are optimal. 

Indeed, as mentioned before, it is easy to check that for the previous range of 
dimensions, u* as in (11.191) is a singular W l,p semi-stable solution — independently 
of the more precise fact of actually being the extremal solution. This remark also 
applies to the following example. 
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Example 1.5. Consider the power nonlinearity f(u) = (l+u) m , with m > p — 1, 
which is studied in [7] and [16] for general domains. Let us define the following 
critical exponent 



( (p-l)n-2y/(p-l)(n-l) + 2-p 4p 

if n > p + 



m cs (p) : = < 



n-(p + 2)-2^g 
+oo if n < p + 



p — 1 

4p 
p — 1 



The results of [7J [16] show that if m < m cs (p) then the extremal solution u* of 
(1.12 Ap ) is bounded. On the other hand, if Q = B 1 and m > m cs {p) then 



u*(r) = r™-^- 1 ) — 1 and A* = I I n 



p \ I mp 



m — (p — 1) J \ m — (p — 1) 

We see that u* e L q (Bi) if and only if 1 < g < n(m — (p — !))/??• Now, if we let 
m = m cs (n,p) then we find that n(m cs (n,p) — (p — l))/p = go and this proves 
the sharpness of Theorem 11.2( c). 

Observe that in the same case m = m cs (n,p) we have 



u (r) = r p v Vh / — 1, 



and, as in the case p — 2, this differs from the pointwise bound (I1.7P for the 
factor |logr| 1//p . Recently, Villegas [29] has proved that, for p = 2, the factor 
| log r I 1 / 2 in (11.71) can be removed. See also [28] for a related improvement of a 
radial Liouville theorem of [1]. 

Remark 1.6. The critical dimension for the boundedness of semi-stable solutions 
determined by n = p + 4p/ (p — 1) tends to +oo as p — > 1 and as p — > +oo. The 
smallest critical dimension corresponds to p = 3, for which n = 9. In the cases 
p = 2 and p = 5 we find the value n = 10. Observe that case (b) in Theorem 11.31 
n = p + 4p/ (p — 1), applies only for those p > 1 such that p + 4p/ (p — 1) is an 
integer. Note also that the critical dimension for the boundedness of all semi- 
stable solutions is bigger than the corresponding dimension for general solutions, 
which is n < p. Therefore, semi-stable solutions enjoy more regularity than 
general solutions. 

Remark 1.7. The following two comments concern the notion of semi-stability 
of solutions. See [7J for a more general nonradial setting. 

(i) If u is a radially decreasing local minimizer as defined after (11. 2p . we claim 
that u is a semi-stable solution. Therefore, Theorem 11.21 applies to every radially 
decreasing local minimizer u G W 1 ' p {Bi). To check this, note that the second 
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variation of the energy functional (jl.2p in Bi at a radially decreasing solution u 
is given by 



p-2 



(P-2) 



Vm 
|Vd 



+ |V£| 2 -</(u)e 2 Ux, 



1.20) 



for every perturbation (gC 1 with compact support in B\ \ {0} (not necessarily 
radially symmetric). It follows that, if u is a radial local minimizer, then (ll.20p is 
nonnegative for every radial £ G C\(B\ \ {0}). Therefore, u is semi-stable in the 
sense of Definition 11.11 since for radial perturbations £, (11.201) reduces to fl 1 . 3 [) . 

(ii) If u is a radially decreasing semi-stable solution of (jl.lj) in the sense of 
Definition ll.l| we claim that (ll.20p is nonnegative for every nonradial £ G C l 
with compact support in Bi \ {0}. That is, the second variation of energy at 
u is nonnegative not only for radial perturbations but also for nonradial ones. 
Indeed, for x G M n we set x = r8, where r = \x\ > 0, and 9 = x/r G dBi = {8 G 
W n : \9\ = 1}. Let £ G C\[B\ \ {0}) (not necessarily radial). We consider the 
spherical averages of £ 2 and define the following radial function: 



\9Bx | 



e 2 (H 



d8 



e 2 (r6 



9B, 



Differentiating the last expression with respect to r and using the Cauchy-Schwarz 
inequality, we find 



dB! 



yV£(r8)] d6= j \UrO)\ 2 d6. 

dB! 



Finally, after some straightforward computations, and using the semi-stability of 
u in the sense of Definition 11.11 for the radial test function (f = tp(r) and that 
C 2 < |V£| 2 in fTOOj) . the claim follows. 

The paper is organized as follows. In Section [2] we prove the main estimates 
on semi-stable stable solutions needed in the rest of the paper. In Section [3] we 
prove Theorem ll.2| while Section H] is devoted to prove Theorem 11.31 



2 Estimates for semi-stable solutions 

This section is devoted to show Lemma |2T3| which contains the key estimate used 
in the proof of Theorem 11.21 

We begin with the following remark on radial functions and radially decreasing 
solutions of (11.11) . 
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Remark 2.1. Given 5 G (0,1), every radial function in W l ' p {B{) also belongs 
(as a function of r = \x\) to the Sobolev space W 1)P {5, 1) in one dimension. As 
a consequence, by the Sobolev embedding in one dimension, u is a continuous 
function of r G [5, 1] , 

\u(l)\ < C n!P \\u\\wi,p(B-L) and Ipll^^^) < C nt p t s\\u\\ w i,p( Bl ) (2.1) 

for some constant C n ^ (respectively, C niPt g) depending only on n and p (respec- 
tively, on n, p, and 5). 

Now, let u G W /1 ' p (i?i) be a radial solution to (11 .ip such that u r {r) < for all 
r G (0, 1). By (12. ip . the right-hand side of (II. ip is bounded away from the origin. 
Hence, by standard regularity theory for the p-Laplacian, u G C^{Bi \ {0}) for 
some (3 G (0, 1). In particular, u T is a continuous function in Bi \ {0} which does 
not vanish and thus (11.11) is uniformly elliptic in compact sets of B\ \ {0}. It 
follows that u G C^(Bi \ {0}). 

Since u r < in B\ \ {0}, we can write (II. ip in radial coordinates as 

- r x ' n d r (r n ~ V| P 'V) = g{u) for r G (0, 1), (2.2) 
and also in the form 

Tii — 1 

- (p - l)\u r \ p - 2 d r u r \u r \ p - 2 u r = g(u) for r G (0, 1). (2.3) 

r 

In order to show our main estimate (stated in Lemma 12. 3j) we need a prelimi- 
nary result on the form of the second variation of the energy for radial solutions. 
This lemma was inspired in the proof of Simons theorem on the nonexistence of 
singular minimal cones in M n for n < 7 (see [U [6] for further details). 

Lemma 2.2. Letu G W 1 ' p {Bi) be a radial solution m5i\{0} of f 1 1.11) satisfying 
u r (r) < for all r G (0, 1) and Q be the quadratic form defined in (11.31) . Then, 

Q{u r r))=J^ K| p j(p-l)|77 r | 2 -^^77 2 j dx, (2.4) 

for every radial function r\ G C\{B\ \ {0}) ; that is, with compact support in 
B 1 \{0}. 

Note that expression ( 12. 4h for the quadratic form Q does not contain any 
reference to the nonlinearity g. This is the reason why our estimates do not 
depend on the nonlinearity g. 

Proof of LemmafZM By Remark EH u G Cf G f {B 1 \ {0}) for some (3 G (0, 1) and 
u satisfies equation ( 12.31) . Let us set H{u) — (p— l)|-u r | p ~ 2 . Let r\ G C\{B\ \ {0}) 
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be a radial function with compact support in B\ \ {0} and c G C X {B\ \ {0}) be 
a radial function. Take £ = cr] G Cl{B\ \ {0}) in (11. 3p to obtain 

Q( cr? ) = / {iJ(n)|V(cr/)| 2 -^(?i)cV}rfa; 

= / {H(u)(c 2 \Vr]\ 2 + Vrf -cVc + r, 2 \Vc\ 2 ) - g\u)c 2 n 2 } dx 

= [ {H{u)c 2 \Vr]\ 2 + H(u)V(r] 2 c) ■ Vc - g'{u)c 2 r} 2 } dx. (2.5) 

Next, we multiply ( 12.3P by 9 r (?7 2 M r r' 1 ~ 1 ), we integrate in r from to 1 and 
use integration by parts to obtain 

= jd r (rfu r r n - l )H(u)d r u r dr - jd r j^-^K| p -V + g{u) | rfu r r n - l dr. 

Using d r (\u r \ p ~ 2 u r ) = H{u)u rr and r n ~ x dr = {dBil^dx, we deduce 

0=1 H(u)d r (r] 2 u r )d r u r dx — [ I g' (u)i] 2 u 2 . — \u r \ p ~ 2 u 2 rj 2 1 dx. 



Therefore 



J {H{u)V{rfu r )-Vu r - g\u)rfu 2 }dx = - J ^—^r] 2 \u r \ p dx. (2.6) 

Taking c = u r in (I2.5p and using (12.61) . we obtain (12. 4L □ 

Now, we use Lemma 12.21 and the semi-stability assumption to establish the 
following result. It is an L p estimate for u r r~ 2a / p in Bi, for certain positive 
exponent a depending on n and p, in terms of the W ' p norm of u. As said 
before, this is the key estimate in the proof of Theorem 11.21 Here as in the 
rest of this section, we assume that n > p. Note that when n < p, we have 
W 1,p (Bi) C L°°(-Bi) and hence solutions are bounded. 

Lemma 2.3. Assume n > p. Let u G W 1)P {B{) be a semi-stable radial solution 
in Bi \ {0} of (II. ip satisfying u r {r) < for r G (0, 1). Let a satisfy 



n — 1 , 

l<a<l + W -■ 2.7 

p — 1 



Then 

a 



\u, 



p r~ 2a dr < — £ II V?;ll p (2 t 



where C ntP is a constant depending only on n and p. 
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Proof. By the semi-stability of u and Lemma 12.21 applied with r\ replaced by rrj, 
we have that 

(n-1) / \u r \ p T] 2 dx < (p- 1) / H P |V(r7/)| 2 cfe (2.9) 

holds for every radial rj e C^(Bi \ {0}). Since 77 vanishes in a neighborhood of 
the origin and u 6 C 1 away from the origin (see Remark 12. ip . we deduce by an 
approximation argument that (I2.9P also holds for every radial Lipschitz function 
rj vanishing on OBi and also in a neighborhood of the origin. 

We now prove that (12. 9p also holds for every radial Lipschitz function van- 
ishing on dBi — but now not necessarily vanishing around 0. To see this, take 
C G C x (SL n ) such that < C < 1, C = in B u and C = 1 in W 1 \ B 2 . Let 
Cj(-) := C(-/S) for every 5 > 0. Replacing rj by rjQ (which is Lipschitz and 
vanishes in a neighborhood of {0}) in ( 12.91) . we obtain 

(n-1) I \u r \ p C, 2 rf dx < (p- 1) I \u T \ p \V{rr]Q\ 2 dx. (2.10) 

JBi JB 1 

Now, using that r\ and V77 are bounded and denoting Ag^s = B 2 s \ B$ we find 
\u r \ p \V{r7]Q\ 2 dx 

= j B \u r \ p {\V(r V )\ 2 C 2 + r 2 V 2 \VCf + CyC s -^(rV)}dx 

< j B \u r \ p \V(r V )\ 2 ( 2 dx + C jf lurrNl^N + ^lCjlVC^lJdo; 

< / |u r | p |V(rr/)| 2 C i 2 rfx + C / |w r | p cix, (2.11) 

where C denotes different positive constants. Since u G W 1,p (Bi), the last term 
in (12.111) tends to zero as 5 goes to zero. Therefore using that Q tend to 1 a.e. in 
Bi as 5 — > 0, (12. lip , and (I2.10p . we obtain by monotone convergence that (12.91) 
holds for every radial Lipschitz function rj vanishing on dB\. 
For a satisfying ( 12. 7ft and e e (0, 1), let 



e _a - 1 for < r < £ 



a Lipschitz function vanishing on dB\. From inequality (12. 9p applied with r\ = r\ e 
we obtain 

(n-1) [ \u r \ p (r- a -l) 2 dx+(n-l)(e~ a -I) 2 [ \u r \ p dx 

i_Bi\B e JB e 

<(p-l)[ \u r \ p {(l-a)r- a -l} 2 dx + (p-l)(e- a -l) 2 [ \u r \ p dx. 

JBi\Be JB e 
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Since n > p it follows that 

(n-l) / H P (r~ a - l) 2 rix < (p-1) / |M r | p {(l-a)r^ Q -l} 2 rfx. 

Jbx\B £ Jb!\B £ 

Developing the squares, using n > p and (12. 7p . we find the estimate 

/ \u r \Pr- 2a dx < ? n '\, 2( 77 / \u r \ p r- a dx. (2.12) 

Jb 1 \b c {n-l) -{a- \) 2 {p - 1) J BABe 

Throughout the proof C n , P (respectively C n , P , a ) denote different positive con- 
stants depending only on n and p (respectively on n, p, and a). 
Now, choose a positive constant C niP>a such that 



(n-l)-(«-l) 2 (p-l) ~2 



r" a < -r- 2a + C n ^ a for all r e (0, f 



Combining this inequality and (12.121) . we obtain 

d/Cc < ^~- xi p o. j I I doc . 

JBiXBe ' ' J Bi\B e 

Now, let e — > to conclude 



it 



r \ p r dx < Cn,p, all ^ M llz / p(ij 1 )- (2-13) 



Bi 



In order to find a more precise expression on how the previous constant C nyPya 
depends on a, we apply (12.131) with the especial choice a = a given by 



1 1 n-l 
ao= 2 + 2V^T G 



, n — 1 



We deduce 

f _ 1 / n-l 

/ K| p r V P -idx<C niP ||V«||^ (Bir (2.14) 
Jbi 

Finally, since r" a < r"( 1+v/(n_1)/(p_1) ) in Bi, (T2TT2|) and (l2~TD lead to the desired 
estimate (I2.8P after letting e — > 0. □ 

3 Pointwise, L g , and VI/ 1 ' 9 estimates 

In this section we prove Theorem 11.21 
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Proof of Theorem \1.2[ By ( 12. ip . \u(l)\ < C ntP \\u\\ w i, P ( Bl y Hence, in view of the 
estimates that we need to prove, it suffices to establish them with u replaced by 
u — u(l) — a positive solution vanishing on dB\. Thus, through the proof, we 
assume 

u>0 = u(l) in Si. 

In case n < p, since u G W 1,p {Bi) ) the Sobolev embedding leads to u G 
L°°(Si) and ||m||z«>(Bi) < C n , p \\u\\ w i, P ( Bl y 

In case n > p, let a satisfy ( 12. 7p . For < t < 1 we have 

U (t) = J _ Mrr -(2«-n+l)/P r (2a-n+l)/ Prfr 

< C n>p (^J \u r \ p r- 2a dxy ^ r p '( 2a - n+1)/p d r y , (3.1) 

by Holder inequality and where p' = p/(p — 1). Using Lemma [2731 (which requires 
n > p) we deduce 

U(t) < C n, P h\\w^(B l} ( f 1 rP '(2a-n + l)/ Pd \ 7 (33) 

{(n-l)-(a-l) 2 (j9-l)}AA / 

for all < t < 1. 

(a) Assume n < p + 4p/ (p — 1). By the remark made above, we may assume 
that n > p. Note that the integral in (13.21) is finite with t = whenever p'(2a — 
n + l)/p > —1, or equivalently 



f 1 rP i(2 a -n + l)/ Pdr < +QO if !^ 



< a. (3.3) 



Recalling (II. 5p and since n < p + 4p/ (p — 1), we can choose a (depending only 
on n and p) satisfying 

n — p In — 1 

< a < 1 1 



2 Vp-1 

In addition, we may take a > 1 — as required in (12.71) . Now, the desired L c 
estimate follows from (13.21) and (I3.3p . 

(b) Assume n = p + 4p/ (p — 1). Let e G (0, 1) and 



In — 1 2p 

a = 1 + V 7 ~ £ = 7 ~ e - 

\ p — 1 p — 1 

Since a satisfies (12 .7\) in Lemma [2731 (13. 2p yields 

u(t) < ^IM^gl) f r -i-^'/P dr 



< gkNkMgO t -2e/p > (3 . 4) 
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for all < t < 1 and < e < 1, where C p = C p , n is a constant depending only in p 
(since here n is a function of p). It follows that u G L q {B{) for every 1 < q < oo. 

In order to prove the pointwise estimate (11.61) . we optimize the right-hand side 
of (13.41) with respect to e by choosing e = log 2| logt| _1 . Note that e is admissible 
since it belongs to (0, 1) if < t < 1/2. With this choice of e, (13. 4ft yields 

u(t) < C p ||m|| W 'i,p(_b 1 )| logt| for < i < 1/2. 

Using this and that u is positive and decreasing, the desired logarithmic estimate 
(USD follows. 

(c) Assume n > p + 4p/ (p — 1) and 1 < q < qo, for q defined as in (11.41) . For 
e E (0,1), let 



n-L 

a = 1 + \l - e. 

p — 1 



By ( 13. 2ft . we have 




1 



1 -I(n-2 A /^-4-p-2+2 e 



< C n , p IMIwn^BO-r* p V Vf-i 7, (3.5) 



£ p 



for < t < 1 , where we have used that — (n — 2 y^^— j- — p — 2 J < since 

n > p + 4p/ (p — 1) — see (ll.5p . 
Now, from (13.51) we obtain 



u q dx< ' p T ( m vp-i )n n - x dt. 



<± 



If we set 

< go 



n-2^^i-p-2 + 3e 

and e > is small enough, the second integral of the previous inequality is finite. 
Hence, u G L (? (S 1 ) for every 1 < q < q . 

Finally, to prove the pointwise estimate ( 11.71) we consider ( 13. 51) and proceed 
as in part (b). Now, we need to make t~ 2£ ' p /e l l p small for given t. We take 
e = log 2 1 log i | — 1 , which belongs to (0, 1) if < t < 1/2. With this choice of e, 
(13.51) leads to 



u(t) < c n , p \\u\\ W i, P{Bl) t ^v n Vp-i p 2 ;|io g t|? 

for t G (0, 1/2). Since u is positive and decreasing, this leads to estimate (11.71) in 
all Bx. 
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(d) Assume g > 0. We prove part (dl). Note that 

r (r n ~ 1 |u P | I '~ 1 ) = -^(r^K^-V) = r^giu) > in B x (3.6) 

and hence r n-1 |w r | p ~ 1 is a nonnegative and nondecreasing function of r. In par- 
ticular 

||r n - 1 k| p - 1 |U»(B 1 ) < \ur(l)r l - (3.7) 

By (13.61) and since u is bounded in S 1 \S 1 / 2 , r n_1 |'u r .| p ~ 1 is a W 1,1 (l/2, 1) function 
of r (indeed W 1,OQ ). By the Sobolev embedding in one dimension, r" --1 ^^" 1 is 
a continuous function up to r = 1. Thus |w r (l)| in (13.71) is well defined and 
bounded. To control it, multiply equation (11.11) by (f £ (r) := min{l, £ -1 (l — r)}, 
where e G (0, 1). This yields 

- / \u r \ p ~ l dx = / g(u)f £ dx. (3.8) 

£ JBi\Bi_ E JBi 

To fully justify (13. 8p . we need to multiply (II. ip by (p E (s, where Q vanishes around 
as in the proof of Lemma [2.31 Then we let 5 — > and use that u G W 1,P (B\) if 
n > p, and that u G C 1 (i?i) and u r (0) = if n < p. 

Now, letting e — > in (13.80 we obtain |ii r (l)| p_1 < C niP ||^(w)||i,i(B 1 ). Thus, by 
(13. 7p we deduce that 

IK-'kMU^Bo < OJbMIUw- (3.9) 

To control || VwUlp^), assume first that n < p. Then, 

r n ~ \u r \ p dr = / (|u r | p r n l ) P -i r p^dr 



o ./o 



1 n-l 



J 

Since the last integral is finite, this and (I3.9P lead to (11.81) . 

In case n > p, we use Lemma 12.31 We take a satisfying (12 .7p and depending 
only on n and p. The lemma gives that 

r- 2a \u r \ p dx (3.10) 

Bi 

< C niP I \u r \ p dx = Cn )P I \u r \ p dx + Cn )P I \u r \ p dx; (3.11) 

here we choose r G (0,1) satisfying Tq 2 " > 2C rhP . Since C n , p j B \u r \ p dx < 
(1/2) j B r~ 2a \u r \ p dx < (1/2) J B r~ 2a \u r \ p dx, we can absorb the first term in 
the right-hand side of (13. lip into (13.101) . and deduce that 

/ \u r \ p dx < [ r- 2a \u r \ p dx < 2C n>p [ _ \u r \ p dx. (3.12) 

JB! JB\ ' J Bi\B rn 
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Note that r depends only on n and p, and thus, since u is decreasing, M(r ) p 1 < 
Cn,p\\u p ^ l \\L 1 {Br Q )- Using this and (13.91) . we have 

I _ \u r \ p dx = C n [ \u r \ p r n - l dr 

JB^Bro Jr 

< C , n ||r n ~ 1 |Mr| p ~ 1 |U°°(Bi) / -u r dr < C ntP \\g(u)\\ LHBl )\\u p ~ l \\^ Bi) . 

Jro 

After using Young's inequality, this bound in B x \ B ro and estimate (13.121) yield 
the desired bound (11.81) in all B\. 

Next, we prove (d2). We consider first the case n < p + 4p/ (p — 1). Observe 
that by part (a) we have that u G L°°(Bi). Hence, applying the regularity theory 
for the p-Laplacian (see [2T]) it follows that u G C 1,/3 (i?i) for some < (3 < 1. 
Moreover, by the L°° estimate of part (a) above and estimate fll.8p . there exists a 
constant depending on n, p, \\u\\ L i( Bl ), and upper bounds in g such that estimate 
(11. 9p holds. This concludes the proof for n < p + 4p/ (p — 1). 

Assume n > p + Ap/{p — 1). We establish parts (d2) and (d3) at the same 
time. Observe that it is enough to prove our estimates in Biu. Indeed, since u 
is decreasing, we have bounds for the supremum of u in Si \ By^ in terms of 
IMIli^). Thus by |2T] we have that u G W l,q (Bi\Bi/ 5 ) for all q < oo. Moreover, 
IMIw^cbiVBi/s) — where C is a constant depending on n, p, \\u\\l^(Bi), and 
upper bounds in g. 

We choose p G (1/4, 1/2) such that 

< -u r (p) = M(1/4) x ~ 4 M(1/2) < 4||u|U<» (1/ 4,i) < C n Ju\\ w i, P{Bl) , (3.13) 

where in the last inequality we have used (12. ip . On the other hand, from (12.21) 
we have 

d r (\u r \ p ~ 2 u r ) = \u r \ p ~ 2 u r — g(u) < \u r \ p ~ 2 u r , (3-14) 

for r G (0, 1). We integrate (I3.14p with respect to r, from t to p, and use (I3.13p . 
to obtain 

f p \u \' p ~ 2 

- \u r {t)\ p - 2 u r {t) < -\u r {p)\ p - 2 u r {p) - (n- 1) / — u r dr 



-J f 1/2 Wrl^ 1 

< C n, P \\u\\ P w l P(Bi) + (n - I) J -y dr, (3.15) 
for all < t < 1/4. Next, we estimate the last integral in ( 13.151) using Holder 
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inequality, to find 



1/2 L. IP-1 /-1/2 

^11 dr = \ u \P-l r -Vc*-n+lW r V*-n+l- P ')/p' dr 



<C n [ I \u r \ p r- 2a dx ) ( r p{2a - n+l - p ' )/p ' dr ) . (3.16) 



Sl/2 / W* 



Let e G (0, 1) and 



/n - 1 

a = 1 + W -£. 

p — 1 



Applying Lemma 12.31 in (13.161) and using that p(2a — n + 1 — p') + p' < since 
n > p + 4p/ (p — 1), we deduce from (13.151) and (13.161) the following estimate 



Wr(*) < 4g IklkMBO t P ^ K (3.17) 



for < t < 1/4. 

Now we use (13.171) to obtain, for q > 1, 



/ |« r | 9 da;<^(||«|| wl , P(Bl) J r P V (3.18) 



If we set 

n-2^-2 + 3e 

the second integral in (13.181) is finite for every e > 0. Hence, w r G L <? (5 1 / 4 ) for 
every q < q±. The estimate (11.91) follows, using in addition estimate (II. 8p . 

Finally, the pointwise estimate (11.101) follows from (I3.17P by choosing e = 
log 4 1 logtl" 1 for < t < 1/4. □ 



4 The extremal solution 

This section is devoted to prove Theorem 11.31 By estimate (11.81) of Theorem 11.21 
now our task is to bound vf^ 1 and f(u\) in L 1 (_B 1 ) uniformly in A G (0, A*). We 
prove this fact using the growth condition (11.131) on / together with the radially 
decreasing property of the minimal solutions U\. These bounds lead to a control 
of u\ in W 1,p (Bi) uniformly in A, by estimate (11.81) . 

Proof of Theorem \1.3[ We start noting that, for A G (0, A*), the minimal solution 
U\ G L°°(Bi) is radially decreasing. This follows from a general result of [5] on 
radially decreasing symmetry, which only requires / to be positive in (0, oo) and 
locally Lipschitz in [0, oo) — see Corollary 1.1 in [9]. We can give, however, a 
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simple proof in our situation of minimal solutions, as follows. First, the minimal 
solution u\ can be constructed by monotone iteration (see [7j), starting with 
u° = 0, by solving —A p u k = A/(ti (r)). Since the right-hand side is given and 
radial, the unique solution u k must be radial. Thus, the limit u\, as k — > oo, is 
radial. Now, integrating our equation —d r {r n ~ 1 \d r u\\ p ~ 2 d r u\) = r n ~ l Xf{u\) > 
in r, from to t G (0, 1), we get that d r u\{t) < 0. 
For A G (0, A*), let p x G (1/2, 1) be such that 

= (2t< A (l/2)r 1 < C^Hu'-'llytB.fl), 

where we have used that u\ is radially decreasing. Since r n ~ 1 \d r u\\ p ~ 1 is increas- 
ing in r G (0, 1) by (13.61) . the previous estimate yields 

iir"- i |9 rMA r- i n i o 0(Bi/2) < pi- i \d r u x { Px )r i < ojk^iili^). (4.1) 

Multiply equation (1.12a iP ) by tp{r) = min{l, 4(1/2 — r) + }, to obtain 
\\Xf(ux)\W(B 1/4 ) < C n , p f' 2 r^dru^dr 

Jl/4 

< C ntP \\u p x ||li(b 1/2 ), (4.2) 

where we have used (14.11) . Note that here we can test the equation with ip, a 
function which does not vanish around the origin, since we know that u\ is an 
energy solution in all of B\. 

Next, we use assumption fll . 13j) on / to ensure, given any 5 > 0, that \f(t) > 
|t p_1 - C s for all t > and A G (A*/2, A*), where C s does not depend on A. This 
combined with (14. 2p leads to 

IK' 1 \\lhb 1/4 ) < C n , P 5\\ uV \ l \W(B 1/2 ) + C s . (4.3) 
Now, since u\ is decreasing in r, we have 

hTXHB^B^) < c n>P ur\i/4) < c n)P \\uT l \W {Bl/A) . (4.4) 

This combined with (14.31) yields, after taking 5 small enough, ||it^~ ||ii(B 1/4 ) < C 
for a constant C independent of A. This bound and the argument used in (14.41) 
done now on B\ \ B\/± leads to a uniform in A bound for llw^ -1 ^ 1 ^)- 

The previous bound gives also a control for ||/(wa)||li(b 1/4 ) ) by (14.21) . Since 
/ is increasing, f(u\) is decreasing in r and the argument used above allows to 
control || /(wa) || uniformly in A. Thus, 

IK^IIncsi) + II/(ua)||h(jji) < c 
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for some constant C independent of A G (0, A*). By estimate (I1.8P of Theorem II. 2 [ 
we deduce a bound for || Vux\\lp(Bi) independent of A. Since u\\ dBi = 0, this yields 
a bound for ||ma||h/ 1 .p(b 1 )- 

Thus, we have that there exists v G Wq' p (Bi) and some subsequence still 
denoted by u\ such that, as A — > A*, u\ — v weakly in W 1,p (Bi), u\ — > v 
strongly in L p (Bi), and u\ — > v a.e. in B\. Since u\ is increasing in A, every 
sequence U\{r) tends as A — > A* to u*(r), by the definition (11.141) of u*. Hence, 
v = u * G Wq' v (B x ). 

Next, we show that u* is an energy solution of (1.12x* lP ). We want to pass to 
the limit in 

/ \Vu x \ p - 2 Vu x -V^dx = [ \f(u x )<pdx for ip G C\{Bx). (4.5) 

JBi JBx 

Since f(u\) increases to f(u*) and 1 1 (iaa) 1 1 x, 1 ( .Bx) is uniformly bounded in A, the 
monotone convergence theorem gives that f(u*) G L l (B{) and that the limit as 
A — > A* of the right-hand side in (14. 5 j) is 

A* f f{u*)ipdx. 

To pass to the limit in the left-hand side of (j4.5j) . note first that we may assume 
that ip is radial — simply by integrating first with respect to the "angles" 9 G dB\ 
in (14.51) . Now, the left-hand side of (14. 5p . up to a multiplicative constant, reads 

/ r n - l \d r u x \ p - l {-d r <p)(r)dr. (4.6) 
Jo 

Using ([3Jl and that the L^O, 1) norm of r n 1 f{u x (r)) is bounded uniformly 
in A, we deduce that r n ~ 1 |<9 r , UA| p ~ 1 is bounded in W rl,1 (0, 1) — the Sobolev space 
in one dimension — uniformly in A G (0,A*). Since this space is compactly em- 
bedded in L 1 (0, 1), we conclude that r n_1 |9 r UA| p ^ 1 converges strongly in L 1 (0, 1) 
to r n_1 |<9 r u*| p ~ 1 . Hence, we can pass to the limit in (14.61) and conclude 

f |W| p - 2 W- Vipdx = X* [ f(u*)tpdx, 
Jb x Jb ± 

for every (f G C^(Bi), proving the claim. 

Clearly, u* is a radially symmetric solution and it is nonincreasing. Hence, 
by regularity, u* G C l (Bi \ {0}) since u* G L°°(Bi \ {0}). Therefore, using that 
u* G W 1,p (Bi) and / > 0, the argument at the end of the first paragraph of this 
proof — or Hopf Lemma (see [21] or Lemma A. 3 in [25]) — , lead to u*(t) < for 
all t G (0, 1). That is, u* is radially decreasing. 

At this point we are able to prove that u* is semi-stable in the sense of 
Definition 11.11 That is, 

[ \*f'(u*)edx< [ ( P -lMr 2 \£ r \ 2 dx (4.7) 

JBx JBx 



20 



for every radially symmetric function £ G C*(i?i\{0}). This follows by passing to 
the limit as A — > A* in the corresponding semi-stability property for the minimal 
solutions ma- On the left-hand sides, we simply, use Fatou's lemma (recall that 
/' > 0). On the right-hand side, we use that £ has compact support in B\ \ {0}, 
and that in such compact set we have C 2,f3 estimates for u\ uniformly in A G 
(0, A*). This holds since |<9 r ii A | > in such compact set, and the same is true for 
|<9 r iiA| uniformly in A G (0, A*), and thus we deal with uniformly elliptic equations. 
Thus \d r u\\ converges uniformly to and the same is true for the quantities 
\d r u\\ p ~ 2 , since \u*\ > — note that the exponent p — 2 could be negative. 

Finally, the regularity statements in the theorem follow as a consequence of 
Theorem OJ □ 

Acknowledgments 

The authors were supported by the Minister io de Education y Ciencia (Spain), 
grant MTM2005-07660-C02-01. The authors are partially sponsored by the Eu- 
ropean Science Foundation (ESF) PESC Programme "Global". 

The second author is also supported by the E.U. RTN program "MULTIMAT" 
MRTN-CT-2004-505226. 

The third author was also supported by CMUC/FCT (Coimbra, Portugal). 

References 

[1] H. Brezis, J.L. Vazquez, Blow-up solutions of some nonlinear elliptic prob- 
lems, Rev. Mat. Univ. Compl. Madrid 10 (1997), 443-469. 

[2] X. Cabre, Boundedness of minimizers of semilinear elliptic problems up to 
dimension four, in preparation. 

[3] X. Cabre, Extremal solutions and instantaneous complete blow-up for elliptic 
and parabolic problems, Contemporary Mathematics 446, American Math. 
Soc. 2007, in: Perspectives in Nonlinear Partial Differential Equations: In 
honor of Hai'm Brezis, pp. 159-174. 

[4] X. Cabre, A. Capella, On the stability of radial solutions of semilinear elliptic 
equations in all ofW 1 , C. R. Math. Acad. Sci. Paris 338 (2004), 769-774. 

[5] X. Cabre, A. Capella, Regularity of radial minimizers and extremal solutions 
of semilinear elliptic equations, J. Funct. Anal. 238 (2006), 709-733. 

[6] X. Cabre, A. Capella, Regularity of minimizers for three elliptic problems: 
minimal cones, harmonic maps, and semilinear equations, Pure and Applied 
Math Quarterly 3 (2007), 801-825. 



21 



[7] X. Cabre, M. Sanchon, Semi-stable and extremal solutions of reaction equa- 
tions involving the p-Laplacian, Comm. Pure Appl. Anal. 6 (2007), 43-67. 

[8] M.G. Crandall, P.H. Rabinowitz, Some continuation and variational methods 
for positive solutions of nonlinear elliptic eigenvalue problems, Arch. Ration. 
Mech. Anal. 58 (1975), 207-218. 

[9] L. Damascelli, B. Sciunzi, Regularity, monotonicity and symmetry of positive 
solutions of m-Laplace equations, J. Differential Equations 206 (2004), 483- 
515. 

[10] J. Davila, Singular solutions of semi-linear elliptic problems, preprint. 

[11] J. Davila, L. Dupaigne, Perturbing singular solutions of the G elf and problem, 
Commun. Contemp. Math. 9 (2007), 639-680. 

[12] J. Davila, L. Dupaigne, I. Guerra, M. Montenegro, Stable solutions for the 
bilaplacian with exponential nonlinearity, SIAM J. Math. Anal. 39 (2007), 
565-592 (electronic). 

[13] J. Davila, L. Dupaigne, M. Montenegro, The extremal solution of a boundary 
reaction problem, preprint. 

[14] S. Eidelman, Y. Eidelman, On regularity of the extremal solution of the 
Dirichlet problem for some semilinear elliptic equations of the second order, 
Houston J. Math. 31 (2005), 957-960. 

[15] P. Esposito, Compactness of a nonlinear eigenvalue problem with singular 
nonlinearity, preprint. 

[16] A. Ferrero, On the solutions of quasilinear elliptic equations with a poly- 
nomial-type reaction term, Adv. Differential Equations 9 (2004), 1201-1234. 

[17] J. Garcfa-Azorero, I. Peral, On an Emden-Fowler type equation, Nonlinear 
Anal. 18 (1992), 1085-1097. 

[18] J. Garcia- Azorero, I. Peral, J. P. Puel, Quasilinear problems with exponential 
growth in the reaction term, Nonlinear Anal. 22 (1994), 481-498. 

[19] N. Ghoussoub, Y. Guo, On the partial differential equations of electrostatic 
MEMS devices: stationary case, SIAM J. Math. Anal. 38 (2007), 1423-1449. 

[20] D.D. Joseph, T.S. Lundgren, Quasilinear Dirichlet problems driven by posi- 
tive sources, Arch. Ration. Mech. Anal. 49 (1973), 241-269. 

[21] G.M. Lieberman, Boundary regularity for solutions of degenerate elliptic 
equations, Nonlinear Anal. 11 (1988), 1203-1219. 



22 



[22] F. Mignot, J. P. Puel, Sur une classe de problemes non lineaires avec non- 
linearite positive, croissante, convexe, Comm. Partial Differential Equations 
5 (1980), 791-836. 

[23] G. Nedev, Regularity of the extremal solution of semilinear elliptic equations, 
C. R. Acad. Sci. Paris Ser. I Math. 330 (2000), 997-1002. 

[24] I. Peral, Multiplicity of solutions for the p-Laplacian, International Center 
for Theoretical Physics Lecture Notes, Trieste, 1997. 

[25] S. Sakaguchi, Concavity properties of solutions to some degenerate quasi- 
linear elliptic Dirichlet problems, Ann. Scuola Norm. Sup. Pisa CI. Sci. 14 
(1987), 403-421. 

[26] M. Sanchon, Boundedness of the extremal solution of some p-Laplacian prob- 
lems, Nonlinear Anal. 67 (2007), 281-294. 

[27] M. Sanchon, Existence and regularity of the extremal solution of some non- 
linear elliptic problems related to the p-Laplacian, Potential Anal. 27 (2007), 
217-224. 

[28] S. Villegas, Asymptotic behavior of stable radial solutions of semilinear el- 
liptic equations in M, N , to appear in J. Math. Pures et Appliquees. 

[29] S. Villegas, Sharp estimates for semi-stable radial solutions of semilinear 
elliptic equations, preprint. 



23 



